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Abstract 

The three-nucleon ground state and the N-d scattering states are obtained us¬ 
ing variational principles. The wave function of the system is decomposed into 
angular-spin-isospin channels and the corresponding two dimensional spatial 
amplitudes are expanded in a correlated polynomial basis. For the scattering 
states, the complex form of the Kohn variational principle is used to deter¬ 
mine the S-matrix. Special attention is given to the convergence pattern of 
the phase-shift and mixing parameters. The calculations have been performed 
using realistic local NN potentials and three-nucleon forces. Important fea¬ 
tures of the method are anomaly-free solutions and the low dimensionality of 
the matrices involved allowing for the inclusion of a large number of states. 
Very precise and stable numerical results have been obtained. 
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I. INTRODUCTION 


Variational principles can be used to obtain the scattering states of few- 
body systems. Among them, the Kohn variational principle has recently 
received renewed attention. In its usual form, the Kohn-Hulthen method is 
applied to calculate the reactance A'—matrix or its inverse AT -1 . In ref Q it 
was pointed out that the presence of singularities in the method was related 
to a zero eigenvalue in the spectrum of the operator H — E. When the oper¬ 
ator H — E is expanded in a real set of basis functions, the associated matrix 
presents a number of eigenvalues close to zero and this number increases with 
the dimension of the basis. In order to avoid these anomalous singularities the 
Kohn variational principle was implemented for the A-matrix or T-matrix 
using proper boundary conditions. In this case the related quantities are 
complex matrices and so the singularities can be shifted to complex energies. 
Occasionally, singularities can also occur in the complex Kohn method, as was 
revealed in ref. [|j]. But they can be eliminated with a judicious choice of the 
non-linear variational parameters. Moreover, the numerical check of the uni- 
tarity of A can be used to detect an anomalous behavior. Applications of the 
complex Kohn variational method have been made in atomic and molecular 
physics for electron-atom and electron-molecule collisions and in nuclear 
physics in nucleon-nucleon (NN) scattering J|]. The principal benefits of the 
method are anomaly-free solutions and the simplicity in the computation of 
the required matrix elements. 

An important problem in variational calculations is the study of the con¬ 
vergence pattern. When the wave function of the system is expanded on a 
certain basis, in general a very large number of elements are needed to get 
fully converging results. This is particularly so in nuclear physics since the 
potential contains a large repulsion at short distances and is strongly state- 
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dependent. Consequently, many states with different spin, isospin and angu¬ 
lar momentum quantum numbers must be included in a single calculation. 

In ref. || the usual Kohn variational principle for the A'-matrix was used to 
describe N-d scattering below the deuteron breakup. The S-matrix was then 
obtained from the relation 

5 = (1 + iK)(l — iK)~ l . (1) 

However, as was pointed out in ref. |J, it is not equivalent to applying the 
complex form of the principle directly to the A-matrix since the above relation 
holds only for the exact matrices. 

In the present paper the complex form of the Kohn variational principle is 
applied to calculating N-d scattering processes. Using the technique of ref. [Q] 
the three-nucleon wave function is decomposed into channels labeled by spin, 
isospin and angular momenta quantum numbers. The corresponding two- 
dimensional spatial amplitudes are expanded in terms of the Pair Correlated 
Hyperspherical Harmonic (PHH) basis ||. The remaining unknown hyper- 
radial functions are taken as a product of a linear combination of Laguerre 
polynomials and an exponential tail. 

Realistic interactions are considered for investigating the flexibility of the 
basis used to describe the structure the system. The Argone AV14 0 nucleon- 
nucleon (NN) potential is used together with three-nucleon interaction (TNI) 
terms. The point Coulomb potential is included in the two proton interaction. 

In section 2 the three-nucleon bound states are calculated as a first application 
of the polynomial basis. Useful information about the number of basis ele¬ 
ments needed to reach convergence and the corresponding numerical accuracy 
is obtained. Section 3 is devoted to the description of N-d scattering using 
the complex Kohn variational principle. Scattering lengths and phase shift 
and mixing parameters have been calculated. In the N-d scattering, states 
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with total angular momentum and parity J n = 1/2" 1 ",3/2 + , 1/2 - ,3/2“, 5/2 
have been considered at the nucleon incident energy Ejy = 3.0 MeV (just 
below the deuteron breakup). Special attention is given to the convergence 
of the elements of the S'-matrix. Variations with respect to the nonlinear pa¬ 
rameters and the fulfillment of the unitarity condition SS' = I are discussed. 
The conclusions and perspectives of the method are given in the last section. 


II. BOUND STATE CALCULATIONS 


Following ref. [0|] the three-nucleon bound state wave function is written 
as a sum of three Faddeev-like amplitudes 


T = i>(yLi,yi) + ip(x j ,y j ) + , 


( 2 ) 


where Xj,y* are the internal Jacobi coordinates which are defined in terms of 
the particle coordinates as 


Xj = Tj - r k , yi = -j=(rj +r k - 2r;) . 


Each i-anrplitude has total angular momentum J J z and total isospin TT Z . 

Using LS coupling it can be decomposed into channels 

N c 

^(xj, y») = ^2^a(xi,yi)y a (jk,i) 


y a {jk,i ) = {[Yi a {xi)Y La (yi)\ Ac 


\si k s , 


a\S c 




aiTT, > 


where Xi,yi are the moduli of the Jacobi coordinates. Each a-channel is 
labeled by the angular momenta £ a , L a coupled to A a and by the spin (isospin) 
sJ a fta) the pair j, k coupled to the spin (isospin) of the third particle s l a 
(t l a ) to give S a (T). N c is the number of channels taken into account in the 
construction of the wave function and can be increased until convergence is 
reached. The antisymmetrization of the state requires that i a + + E k be 

odd, while the parity of the state is given by l a + L a . 


( 3 ) 


( 4 ) 

( 5 ) 


4 


The two-dimensional spatial amplitude in each channel is expanded in 


terms of the PHH basis: 

£«*(p) ( 2 ) 4 a,La (&)I , 

. K 

where the hyper spherical polynomials are 

( 2 )p^a,ia(^.) = (sin (fii) La (cos ^)j)^ Q P r f“ +1/2 ^ Q+1/2 (cos 2 </>j) . 

Here, N^ x,La is a normalization factor, is a Jacobi polynomial and 

K = £ a + L Q + 2n is the grand orbital quantum number which runs from its 
minimum value Kq = £ a + L a to its maximum selected value K a . Therefore, 
the number of hyperradial functions per channel is M a = (K a — K$)/2 + 1. 
The inclusion of the pair correlation function f a (xi) in the expansion of eq. (] 6 |) , 
as has been discussed in ref. ||, accelerates the convergence taking into ac¬ 
count the correlations introduced by the strong repulsion of the NN potential. 
It is obtained solving a two-body Schrodinger-like equation. 

In the present work the hyperradial functions u c ^(p) are taken as a product 
of a linear combination of Laguerre polynomials and an exponential tail: 

Uk(p) = E4, m £ mWffl t P(-z) > 

m 

where z = 'yp and 7 is a nonlinear variational parameter. Let \a,K,m > 
be a totally antisymmetric element of the expansion basis, where a denotes 
the channel including the angular-spin-isospin dependence and the correla¬ 
tion function, and I\ , m are the indices of the hyperspherical and Laguerre 
polynomials, respectively. In terms of the basis elements the wave function 
(| 2 |) can be writen as 

^ = J2 A K,m\^K, m > . 
a,K,m 

The problem is to determine the linear coefficients AJ <m . To this aim, the 
Hamiltonian considered here includes the AV14 NN potential and the three- 
body Tucson-Melbourne (TM) force The cutoff parameter of the TM 


^a{Xi,Vi) = p la+La f a (Xi) 




potential is chosen to be A = 5.13/U (p is the pion mass) in order to repro¬ 
duce the 3 H binding energy. To evaluate Coulomb (C) effects the AV14+C 
potential model is also analyzed. The choice of the interaction model allows 
for rigorous comparisons with other techniques and it reproduces reasonably 
well the structure of the three-nucleon system. 

The wave function and the energy of the system are obtained by solving 
the following generalized eigenvalue problem 

y A^, m , < a, K, m\H — E\a ', K ', m' >= 0 . 

a' ,K' ,m' 

In the latter equation the dimension of the involved matrices is related to three 
indices: the number of a-channels N c . the number of hyper spherical polyno¬ 
mials for each channel M a and the number of Laguerre polynomials included 
in the expansion of the hyperradial functions (| 8 |). To analyze the convergence 
properties of the expansion, the indices K , m have been increased until stable 
results are reached. Beginning with those channels having minimum grand 
angular quantum number Kq = 0, new channels with higher Kq values are 
successively included. The ground state has total angular momentum and 
parity J n = l/2 + and total isospin T = 1/2. The number of channels with 
Kq < 4 is N c = 18 and with Kq < 6 is N c = 26. In refs. ||11 it was pointed 
out that the number of hyperspherical functions per channel needed when 
using the PHH is small: typically 6 to 8 functions in those channels with the 
deuteron quantum numbers (t a = 0,2;s^ = 1 \t l £ = 0), and 4 to 6 in the 
other ones are sufficient for a four digit accuracy. 

The matrix elements of eq. ( 0 ) have been calculated numerically using a 
Gauss formula for integrals in the variable pt = xi ■ yi, a Chebyshev Lobatto 
formula for integrals in the variable cos 2<pi and a Lagrange formula for the 
integrals in the hyper radius p [9]. 

Special attention has been paid to the study of the convergence with the 


( 10 ) 
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Laguerre polynomials and the nonlinear parameter 7 . The results for the 
three potential models are indicated in fig. 1. In each panel three curves are 
given corresponding to the number of channels N c = 14,18, 26 (the ordering 
of the channels is given, for instance, in ref. |H|). The binding increases with 
the number of channels, but changing N c from 18 to 26 produces an increment 
of about 1 keV. Higher channels are strongly suppressed due to centrifugal 
barrier effects. The N c = 26 results for the binding energy, kinetic energy and 
S-, D- and P —wave occupation percentage are given in table 1. High preci¬ 
sion results from refs. @0 are also shown for the sake of comparison. An 
extremely close agreement has been obtained in all the calculated quantities. 
In the final calculation we have N c = 26, m = 16 and 7 = 2.0 fin -1 . The 
number of basis elements was 1900, which is the dimension of the matrices in 
the generalized eigenvalue problem. 

The variation of the binding energy B with the nonlinear parameter 7 
and with the number of Laguerre polynomials m = 4,8,12,16 is presented in 
fig.2. It can be seen that the dependence on 7 is appreciable only for rather 
low m values. On the other hand, the convergence with m is very fast and 
for m > 12 the results are already accurate at the 0.1% level. The number 
of basis elements is large enough to guarantee a weak dependence with the 
nonlinear parameter. 

The results obtained so far show the great flexibility of the basis for 
reproducing the structure of the system, even when the AV14+TM and 
AV14+Coulomb potential models are used. There are no difficulties in ap¬ 
plying the polynomial basis to other realistic NN interactions such as the 
Argonne AV18 potential, Bonn potential, Nijmegen potential, etc. For those 


interactions expressed in momentum space the matrix elements of eq.(lC) can 
be calculated after performing a Fourier transformation the basis. 
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III. THE COMPLEX KOHN VARIATIONAL PRINCIPLE 


Following ref. |5|, the wave function corresponding to an N-d scattering 
state is written as a sum of two terms 


■ 


( 11 ) 


The d'c; term describes the system when the three-nucleons are close to each 
other. For large interparticle separations and energies below the deuteron 
breakup it goes to zero, whereas for higher energies it must describe three 


outgoing particles [12]. The second term is the solution of the Schrodinger 
equation in the asymptotic nuclear region, where the two incident clusters are 
well separated. It can be written as a linear combination of functions of the 
form 


^Lsjixiiyi) =T,i a =o, 2 Wi a {xi)TZ x L {yi) x 

{ [[*!« (*i)4*]l ® «*] s ® [Va?] TT Z , 


where wi a (xi) is the deuteron component in the state l a = 0,2 and L is the 
relative angular momentum of the deuteron and the incident nucleon. For 
A = R, lZ\{yi) coincides with the regular solution Fj^yi) of the two-body 
N-d Schrodinger equation in absence of nuclear interaction. For A = /, it 
is set equal to Gi(yj), the product of the irregular solution Gl{Ui ) and a 
regularizing factor which goes to zero as yf +1 , when j/j —► 0, and approaches 
unity for large values of y*: 


( 12 ) 


G L (y) = (l-e-t rNd ) L+1 G L (y), 


( 13 ) 


where r^d = (2/\/3 )y is the nucleon-deuteron separation. Hereafter the 
parameter ^ is given the value £ = 0.25 fm _1 |5]. The asymptotic wave 
function can also be constructed in terms of the physical ingoing and outgoing 
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solutions by linearly combining the above functions. As in ref. j|] we can define 
a general asymptotic state 


ttlsAxiiYi) — ^LSj( x i> y*) + Yj J £ LULL'S' j( x j>yi) > (14) 

L' S' 

with the following asymptotic functions 

^LsA^i^i) = «00^LSj( x oy*) + u Olll£sj( x i,yj) , (15) 

Hisj( x *>y*) = uio^LSji^-^yi) + u ntt I Lsj(xi,yi) , (16) 


Four different choices of the matrix C can be obtained in correspondence to 
the following definitions for the matrix u. 

(iA 


C = K for u = 


C = K 1 for u = 


\°V 
foil 


jC = S for u = 


V 10 / 

j -A 

v 1 J 

fio) 


C = —7 tT for u = 


(17) 


(18) 


(19) 


( 20 ) 


V* V 


Accordingly, the asymptotic states have been normalized in order to satisfy 
the Wronskian 


W(yP LS j,Vt l LSJ ) = uoiuio - u 00 un = —det (u) . 


( 21 ) 


The three-nucleon scattering wave function for an incident state with rel¬ 
ative angular momentum L, spin S and total angular momentum J is 


*1sj = Y [^c( x i,yi) + HL5j( x *>y0] , 

i=l,3 

and its complex conjugate is A variational estimate for the matrix C 

can be obtained from the generalized Kohn variational principle 


(22) 
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(23) 


\ 3 L S lv ] = J £ S LL' - ^u)^ + lsj\H ~ E\*+S'j) • 

The trial parameters in the wave function ’b/^j are varied in order to obtain 
a stationary value of the above functional. 

In the present work the calculations have been restricted to energies be¬ 
low the three-body breakup. Therefore, the 'be term goes to zero when the 
hyperradius p —» oo, and it has been expanded in terms of the polynomial 
basis introduced in the previous section. 


*c= E A K,m\ a,K,m>. (24) 

a,K,m 

With this choice of ’be, the variation of the diagonal functionals with respect 
to the linear parameters leads to the following linear system 

E A K',m' < a ) K; m\H — E\a ', K\ m! >= D^ SJ (a, K, m) , (25) 

a' ,K' ,m' 

with the two different inhomogeneous D terms corresponding to A = 0,1 

D x LSJ (a,K,m) = E < a,K,m\H - ^|^isj( x j, Yj) > • (26) 

3 

The first order solution of the matrix L is obtained by solving the following 
algebraic equations 


E J r SS" vS'S" 
l -LL" a -L'L" 

L",S" 


= Y, 


SS' 


LL' > 


with the coefficients X and Y dehned as 


(27) 


Xlu = (Vlsj + *lsj\H - mi'S'j) , ( 28 ) 

Ylv = (n° LS j + *lsj\H ~ E\n° L , SIJ ) , (29) 


where 4 ’^SJ th e solution of the set of eqs.(25) with the corresponding in¬ 
homogeneous term. The second order estimate is obtained from the 

first order solution using eq.( 
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The calculations has been performed applying eq.(^3|) to the 5-matrix. 
The unitarity of the 5-matrix, 55^ = I, has been checked to verify the 
completeness of the basis and to detect eventual spurious solutions. The 
reactance it-matrix, which is real and symmetric, can be obtained from the 
relation 


K = (ttoi + uii£)(it 0 o + u 0 iC) 


-l 


In an N-d process, the total energy of the system is E = —B^ + Eq with 
Bd the binding energy of the deuteron and Eq the incident nucleon energy. 
As a first application of the method, the zero energy process at Eq = 0 is 
studied. In this case only 5-wave scattering has to be considered, with two 
possible states, J = l/2 + and J = 3/2 + . The Hamiltonian of the system 
includes two and three-nucleon forces. As for the bound state, the AV14, 
the AV14+TM and AV14+C potential models are utilized. After fixing the 
number of channels N c , eq.(25) has been solved with increasing values of K, m. 
States up to Kq < 6 have been taken into account, which means N c = 26 for 
J = l/2 + and N c = 44 for J = 3/2“*". The second order results for the doublet 
2 a and quartet 4 a scattering lengths and the difference |55^ — /1 are given in 
table 2 for a few values of the number m of Laguerre polynomials. In this 
process the 5-matrix is a scalar, therefore to check the unitarity condition 
the square of the modulus of a complex number must be calculated. The 
difference between the first and second order results measures the error in the 
construction of the wave function. To give an idea, the relative difference is 
0.1 for m = 4, but it is 10 -4 for rri = 20. The accuracy of the solution is also 
reflected in the unitarity condition which is satisfactorily fulfilled, indicating 
the absence of spurious solutions Q. 

N-d scattering at zero energy has been studied using similar potential 
models in ref. Il3], solving the corresponding Faddeev equations in configu- 


( 30 ) 
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ration space, and in ref. [f|] applying the Kohn variational principle to the 

matrix. The corresponding results are given in the table for the sake of 
comparison. The differences can be attributed to working in a different Hilbert 
space. For example, in ref. [ 131 a truncation of the potential retaining those 
NN partial waves with j < 4 has been considered, whereas in ref. |^] a smaller 
number of channels in the expansion of the wave function has been utilized. 

The calculations reported the table 2 have been done with a value of the 
nonlinear parameter 7 = 1.5 fm . The results for 2 a when 7 is varied around 
this value are shown in the first three panels of fig.3, for the three potential 
models under consideration. A similar behavior is obtained for 4 a, with the 
conclusion that the results have only a slight dependence on 7 when enough 
polynomials are included. It is interesting to note that in this scattering 
calculation the convergence is faster around the value of 7 = 1.5 fm -1 , whereas 
in the bound state case this occurs around 7 = 2.0 fm -1 . The lower value 
of 7 is motivated by the fact that the expansion basis must construct the 
scattering wave function (^) also at intermediate particle separations where 
the term is not yet the correct solution. Finally, in the last panel of 
fig.3 the convergence with the number of channels is displayed. Changing the 
number of channels from N c = 18 to N c = 26 produces a variation of the 
scattering length less than 10~ 3 fm. 

For positive incident energies, S-wave and P —wave scattering have been 
studied. Accordingly, five different states can be constructed with J = 
l/2 + , 3/2 + , 1/2“, 3/2“, 5/2”. When J = 1/2 the corresponding S-matrix 
is a 2 x 2 matrix, in all the other cases it is a 3 x 3 matrix. For each J n state 
all channels with Kq < 6 have been included. It was checked that higher 
channels give negligible contributions. The case considered here corresponds 
to a nucleon incident center of mass energy of E$ = 2.0 MeV (£ 7 ^ = 3.0 
MeV), slightly below the deuteron breakup. The results of main interest are 
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presented in figs.4, 5. In fig.4 the two 5-wave phase shifts and the five P~ 
wave phase shifts are shown for different numbers of Laguerre polynomials 
m, whereas the corresponding mixing parameters are presented in fig.5. The 
solid curve corresponds to the AV14 potential model, the dashed line to the 
AV14+C and the dotted line to the AV14+TM potential models, respectively. 
Stable results are obtained for m > 8 , especially for the mixing parameters. 
An increase of m from 16 to 20 gives no appreciable contribution. More 
interesting is the constant difference between the AV14 and the other two po¬ 
tential models thereby providing a fast method for evaluating Coulomb and 
TNI effects. 

For a more detailed analysis of the convergence properties of the correlated 
basis in conjunction with the complex Kohn principle, the 5-matrix elements 
5££/ are given for the J = 3/2“ state in table 3. This case illustrates the main 
aspects of the convergence with similar patterns in the other four cases. With 
the condition Kq < 6 the number of channels allowed are N c = 36, which 
suffices for a four digit accuracy. In the table, the second order solutions 
for the three phase shifts and the three mixing parameters of the state, are 
given in columns 2-7. The absolute value of the diagonal elements of the 
matrix 55^ — 7 are also shown in the last three columns. The 5-matrix is 
ordered following the convention of Seyler [|b|], therefore the first diagonal 
element corresponds to the A F state, the second one to the 2 P state and the 
third one to the 4 P state. For all the six parameters the convergence is quite 
satisfactory and the unitarity condition is fulfilled. The complicated structure 
of the TNI is reflected in a slightly worse, although satisfactory fulfillment of 
this condition. The behaviour of the phase shift and mixing parameters with 
the nonlinear parameter 7 is similar to that found for the bound state and 
for the scattering lengths. Again a faster convergence is obtained around the 
value 7 = 1.5 fm _1 as in the preceding case. The above analysis can be 
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repeated for states with different J and parity values, yielding in all cases 
quite accurate results. 


IV. CONCLUSIONS 

In the present paper a detailed numerical analysis of the three-nucleon 
bound state and of the N-d scattering have been presented. For the N-d 
process the complex Kohn variational principle has been utilized, limiting 
the study to energies below the deuteron breakup. A correlated polynomial 
basis with exponential tail has been used to expand the internal part of the 
three-nucleon wave function. The Hamiltonian of the system includes realis¬ 
tic two- and three-nucleon interaction terms and the point Coulomb poten¬ 
tial. Explicitly, the Argonne AV14 potential and the Tucson-Melbourne 7 r- 7 r 
three-nucleon force have been considered. Two different motivations have 
been pursued. First of all, the variational method investigated allows for an 
accurate description of the scattering states. Singularity-free solutions have 
been obtained and the unitarity condition for the S'—matrix is satisfied. The 
second motivation was the study of the convergence pattern related to the use 
of the correlated polynomial basis. This basis presents a fast rate of conver¬ 
gence and a good stability with the nonlinear parameter 7 in the exponential 
tail. Binding energies, scattering lengths and phase-shift and mixing param¬ 
eters have been calculated with an error of less than 0.1%. The number of 
basis elements is reasonably small and the related dimension of the matrices 
is kept low. 

The analysis of the convergence was studied paying attention on three 
indices. For a given J 71 " state, all channels corresponding to Kq = £ a + L a < 6 
have been taken into account. At the energies considered here, channels with 
higher values of Kq are strongly suppressed due to centrifugal barrier effects. 
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The number of correlated hyperspherical functions for each channel has been 
increased until a complete stability in the calculated quantities was obtained. 
In this analysis the bound state calculations, presented in section 2, and 
previous experience from ref. || was utilized. The conclusion is that 8 hyper¬ 
spherical functions in those channels with the deuteron quantum numbers and 
4 to 6 in the others suffice for the desired accuracy. Finally, the convergence 
with respect to the number of Laguerre polynomials was carefully studied. 
Increasing m from 16 to 20 does not produce any numerical improvement and 
the convergence is already achieved with m = 12. When the maximum num¬ 
ber of basis elements is used, the relative difference between the first order 
and second order results for the 5-matrix is 10 -4 and the unitarity condition 
is fulfilled up to 10~ 9 (without TNI) and 10 -6 (with TNI). The difference 
reflects the complicated structure of the system when the TNI is present. 

With the method presented here, it will be possible to give accurate the¬ 
oretical predictions for different observables such as the N-d differential cross 
section, vector and tensor analyzing powers, etc. In particular, it is possible to 
use the present method for evaluating TNI effects in polarization observables 
where important differences have been observed between theoretical calcula¬ 


tions and the available experimental data [15]. In conclusion, the extension 
of the method to treat the breakup channel is of great interest. The complex 
Kohn variational principle is well suited for including outgoing boundary con¬ 
ditions in the internal part of the wave function to describe the three outgoing 
nucleons. A first application of the method including the breakup channel is 
given in ref. [O] and is at present subject of intense investigations. 
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Table captions 


Table 1. Binding energy and kinetic energy (in MeV), and S-, D— and 
P —wave occupation probabilities for the three potential models under consid¬ 
eration. 

Table 2. The doublet and quartet scattering lengths (in fm) for different 
values of the number m of Laguerre polynomials and for the three potential 
models considered. The unitarity condition is checked in terms of the absolute 
difference \SS^ — I\. 

Table 3. Phase shift and mixing parameters corresponding to the state 
J = 3/2“ as functions of the number rn of Laguerre polynomials. In the last 
three columns the unitarity of the S-matrix is checked reporting the diagonals 
elements of the matrix difference | SS' — 1 1. 
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Figure captions 


Figure 1. The convergence of the binding energy (in MEY) is shown as a 
function of the number m of Laguerre polynomials. The solid line corresponds 
to the number of channels N c = 26 in the expansion of the wave function; the 
dot-dashed line to N c = 18 and the dashed line to N c = 14. 

Figure 2. The calculated value of the binding energy (in MeV) as a function 
of the nonlinear parameter 7 . The different curves in each panel correspond 
to increasing the number m of Laguerre polynomials. The long dotted line 
corresponds to m = 4, the short dotted line to m = 8 , the dashed line to 
m = 12 and the solid line to m = 16. 

Figure 3. The doublet scattering length (in frn) as a function of the non¬ 
linear parameter 7 (first three panels) and for increasing number of channels 
(fourth panel). In panels 1 — 3 the curves are given for increasing number m 
of Laguerre polynomials, as in fig. 2. In panel 4 the convergence with respect 
to the number of channels is shown for the three potential models considered. 

Figure 4. Convergence of the phase shift parameters (J = 3/2“) with 
respect to the number m of Laguerre polynomials. For the sake of clarity the 
AV14+C is given separately in panels 4 and 6 . 

Figure 5. Convergence of the mixing parameters (J = 3/2~) with respect 
to the number m of Laguerre polynomials. 
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TABLES 


potential 

B(MeV) 

T(MeV) 

Ps(%) 

Pd (%) 

Pp(%) 

AV14 

7.684 

45.678 

90.957 

8.967 

0.076 

ref. 1C] 

7.684 

45.677 

90.956 

8.968 

0.076 

AV14+C 

7.033 

44.813 

90.994 

8.931 

0.075 

ref. 1C] 

7.033 

44.812 

90.993 

8.932 

0.075 

AV14+TM 

ref. 11] 

8.486 

8.486 

49.337 

49.340 

90.577 

9.262 

0.162 


TABLE I. 
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AV14 

m 

2 a (fm) 

|SSt - I\ 


4 a (fm) 

|SSt - I\ 

4 


1.426 

10~ 3 


6.399 

10 -4 

8 


1.203 

10" 4 


6.381 

10 -6 

12 

1.190 

10 -6 


6.379 

10 -7 

16 

1.189 

10~ 8 


6.379 

10“ 8 

20 

1.189 

10~ 9 


6.379 

10“ 9 

ref. § 

1.196 



6.380 


ref. 


1.204 



6.380 



AV14+TM 

m 

2 a (fm) 

|SSt - 1\ 


4 a (fm) 

|SSt - 1\ 

4 


0.8541 

10" 2 


6.391 

10~ 2 

8 


0.6016 

10~ 4 


6.373 

10" 4 

12 

0.5869 

10- 7 


6.371 

10“ 7 

16 

0.5858 

10" 6 


6.371 

10 -6 

20 

0.5857 

10~ 6 


6.371 

10 -6 


AV14+C 

m 

2 a (fm) 

|55t - I\ 


4 a (fm) 

\S& - I\ 

4 


1.369 

10~ 2 


14.002 

10~ 2 

8 


0.9656 

10~ 4 


13.790 

10~ 4 

12 

0.9431 

10~ 5 


13.774 

10 -6 

16 

0.9414 

10' 8 


13.773 

10“ 8 

20 

0.9411 

IQ - 9 


13.773 

10 -9 

ref. § 

0.954 



13.779 


ref. 


0.965 



13.764 



TABLE II. 
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AV14 


m 

^3/2— 

2 ^3/2- 

4 P 3 /2- 

7.3/2- 

e 3/2- 

C3/2- 


1 

H— 

Co 

Co 


4 

0.9403 

-7.488 

26.19 

-3.656 

-2.712 

-0.2715 

10" 4 

10 -3 

10” 2 

8 

0.9425 

-7.216 

26.38 

-3.794 

-2.755 

-0.2592 

10 -6 

10" 7 

10 -5 

12 

0.9427 

-7.198 

26.39 

-3.804 

-2.760 

-0.2580 

10- 7 

10- 7 

10" 6 

16 

0.9428 

-7.197 

26.40 

-3.804 

-2.761 

-0.2578 

10- 7 

10" 7 

10“ 8 

20 

0.9428 

-7.196 

26.40 

-3.804 

-2.761 

-0.2578 

10" 9 

10" 9 

10“ 8 


AV14+TM 

m 

4 ^3/2- 

2 P 3 /2- 

4 P 3/ 2- 

7.3/2- 

e 3/2- 

C3/2- 


1 

H— 

CO 

Co 


4 

0.9408 

-7.456 

26.47 

-3.681 

-2.816 

-0.2588 

10" 4 

10 -3 

10” 2 

8 

0.9430 

-7.183 

26.68 

-3.819 

-2.859 

-0.2465 

10 -6 

10" 6 

10- 4 

12 

0.9432 

-7.166 

26.69 

-3.827 

-2.865 

-0.2448 

10- 7 

10 -6 

10 -6 

16 

0.9433 

-7.163 

26.69 

-3.829 

-2.866 

-0.2448 

10- 7 

10 -6 

10" 6 

20 

0.9433 

-7.163 

26.69 

-3.830 

-2.866 

-0.2448 

10" 8 

10 -6 

10" 6 


AV14+C 

rn 

4 ^3/2- 

2 P 3 /2- 

4 P 3 /2— 

73/2- 

e 3/2- 

C3/2- 


h — 1 

1 

H— 

co 

co 


4 

0.8705 

-7.353 

24.45 

-3.273 

-2.295 

-0.3301 

10- 4 

10 -3 

10 -2 

8 

0.8718 

-7.175 

24.61 

-3.356 

-2.323 

-0.3350 

10" 6 

10" 8 

10" 5 

12 

0.8720 

-7.167 

24.62 

-3.360 

-2.327 

-0.3219 

10" 6 

10- 7 

10" 6 

16 

0.8721 

-7.166 

24.63 

-3.362 

-2.327 

-0.3219 

10" 6 

10- 7 

10" 8 

20 

0.8721 

-7.166 

24.63 

-3.362 

-2.327 

-0.3219 

10" 8 

10" 8 

10" 8 


TABLE III. 
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nonlinear parameter 7 (fm 1 ) 
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1 1.3 1.5 1.8 2 

nonlinear parameter 7 
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